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Abstract 

An external description for aperiodically sampled MIMO linear sys- 
tems has been developed. Emphasis is on the sampling period sequence, 
included among the variables to be handled. The computational proce- 
dure is simple and no use of polynomial matrix theory is required. This 
input/output description is believed to be a basic formulation for its later 
application to the problem of optimal control and/or identification of lin- 
ear dynamical systems. 

Keywords: Balancedness, Bit-string model. Combinational generator. 
Design rules. 

1 Introduction 

There are two different ways of describing dynamical systems: 

(i) by means of input/output relations; 

(ii) by means of state variables. 

In the classical or frequency-domain approach, systems are described by 
transfer functions which reflect just the external or input/output properties of 
the system. However, this mode of description entails some difficulties concern- 
ing stability and realization [I], [2]- 

The modern or time-domain approach turns around the axiomatic concept of 
state. The method is exact in defining the notion of dynamical systems and also 
describes all internal couplings among the system variables [3]; [!]■ Nevertheless, 
the procedure became somewhat disappointing due to the necessity of finding 
state- variable models and to the implicit assumption that all state- variables are 
accessible for direct measurement. This assumption is justified in mechanical 

''Work supported by Ministerio de Educacion y Cicncia (Spain). 
First version published in IEEE Transactions on Automatic Control. Volume 33, No. 4, pp. 
381-384. April 1988. 



or electrical systems but it is not generally satisfied for plants in chemical, gas, 
paper, and other industries. 

These considerations were responsible for the comeback of transfer function 
methods [5]- [7]. 

On the other hand, the enormous increase in the use of digital computers in 
process control has stimulated studies in the field of discrete systems for both 
types of representation. See [8]- [10] and also the above mentioned references. 
All of them are concerned with constant sampling period, which is convenient 
for the simplicity of implementation and mathematical treatment. However, the 
general case of aperiodic sampling is a priori capable of more favorable solutions 
to the problem of control and/or identification of dynamical systems, and it is 
also feasible with modern time-sharing equipment. 

In this work, an input / output modeling technique for aperiodic sampling lin- 
ear systems has been developed. The external description includes the sampling 
sequence among the variables to be handled. The system is described by in- 
put/output data according to the actual experimentation conditions. Although 
the multivariable case is covered, the complexity of the polynomial matrix the- 
ory is avoided. 

The procedure is believed to be a basic formulation for its later application 
to the synthesis of linear control systems sampled in an aperiodic way, since 
most of these techniques for nonperiodically sampled systems rely exclusively 
on the state-space equations [II]-[I3]- 

2 Basic Assumptions 

Our discussion is restricted to the following: 

(a) linear time-invariant multivariable dynamical systems of finite order; 

(b) systems whose transfer function is a p x m matrix (m-inputs, p-outputs), 

where the different entries are strictly proper rational functions. 

We end this preliminary section with the following statement. 
Statement: Let (Gi) be a family of vector functions 

Gi;\R'' -^\R" GieC°°{\K'',\R'') {I ^ 0,1, . . . ,n) 

C°°(|i?", being the set of infinitely differentiable functions on \R". If 
the following conditions are verified: 

(a) there exists an integer r < n such that the elements {Gr{z)) are linearly 
independent for all z G |-R" 

(b) there exists an integer k > r such that {Gr{z)) depends linearly on {Go{z), . . . 

,Gr{z)) 

then, there are functions /o, /i, . . . , /„ G C°°(|i?", \R) such that the following 
expression holds: 
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(1) 



/=0 



The previous result is a direct consequence of the Cramer Rule; for more 
details see [TE\. 



3 External Description for Nonperiodically Sam- 
pled Linear Systems 

3.1 Input/Output Modeling Technique 



Let H{s) be the matrix transfer function of a linear time- invariant multivariable 
system. 

His) = (Hrgis)) (r = 1, . . . ,p), (g = 1, . . . , m) 
let us rewrite H{s) as 

N{s) 



His) 



d(s) 



where 



(2) 
(3) 



dis) 



dis' 



is the least common multiple of the denominators of the entries of H{s). 
In the time domain, the impulse response h[t) can be written as 



h{t) = (hrgit)) = 



hi{t) 



hp{t) 



(4) 



(5) 



where the rth row can also be written in matrix form by means of the triad 
{A,C,Br). In fact, 



hr{t) = C exp{At)Br (r = l,...,p) 
with A — a bottom-companion matrix with last row 

— [dm dn-1, ■ ■ ■ ,di] 



(6) 
(7) 



C= (l,0,...,0)ix„ 



(8) 



Br 



/ KM 

Km 



hrm{0) 



(0) / 



(9) 
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Remark that the column vectors of Br correspond to the n-first Markov 
parameters of the scalar impulse responses hrq{t) (g = 1, . . . , m). 

It should be noticed that the triad {A,C,Br) leads us naturally to the ob- 
servability canonical realization from the vector impulse response hr(t). 

From this triad, we are going to define a family of vector functions Gi; 
{l = 0,...,n) given by 

Gjiz) = C expiA{zi + ... + z,)) (j = l,...,n) (10) 

Go{z) = G (11) 

with 

z= (zi,...,z„) e (12) 

Thus, the vector impulse response hr{t) can be written in terms of these 
functions as 

hr{zi + ... + Zj)^Gj{z)Br (r=l,...,p). (13) 

From an analytic viewpoint, the functions Gi belong to C°°(|i?", ji?") as 
composition of G°° functions. 

It has been proved [T3] that there is an open interval / of |i? such that the 
vectors {Go{z),Gi{z), . . . ,G„_i(2)) defined as before are linearly independent 
for each z€lxlx...xl — I^. 

In this case, it is easy to see that for the new domain /" the conditions a) 
and b) in the previous statement hold. Hence, there will be functions fi(z) G 
C'^il'^, \R) {1 = 0,..., n) such that 

n 

Y,fn-i{z)Gi{z) = -izel^. (14) 

Thus, the functions (/i(z), . . . , fn{z)) can be obtained by solving a compat- 
ible system of linear equations. (For simplicity {fo{z) = —1). 
In fact, the general form of these functions is 

, , , Det{G',{z),...,G'^{z),...,G'^_,iz)) 

fn-i{z) = n ( \ 7^' TY\ (/ = 0, . . . , n - 1) (15) 

-Dei(Go(z),...,G„_i(z)) 

(' denotes the transpose) where the numerator is the determinant obtained from 
the matrix {Gq{z), . . . , G„_]^(z)) by replacing the (Z-l-l)th column by the column 
vector G„(z). 

Now, we multiply both sides of (14) by 

exp{Az*)Br (r = l,...,p) (16) 
with z* taking successively the values 
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I 

z* = -Y.^i (/ = l,...,n) (17) 

i=l 

and we get in each case 

C exp{A{zi+i + . . . + z„))Br 

n—l 

= X! /»(-^)'^ expA{zi + i + . . . + Zn-i))Br 
i=l 

(-1 

+ J2 fn-i{z)C exp{-A{z^+i + ... + Zl))Br (18) 

1=0 

(Z = l,...,n),(r = l,...,p). 

Finally, we define 

9n-i{z) = Y,fn-i{z)C exp{-A{z^+, + ... + Zi))Br (19) 

= {g^^U^),...,g:,-ll{z)). 

At this point, we can identify the components of z with the elements of the 
samphng period sequence. 
In fact, 

Zn-i=tk-i-tk-i+i=Tk-i (Z = 0,...,n-1) (20) 

where tk-i are the sampling instants and Tk-i the length of the sampling 
intervals. 

Thus, at an arbitrary sampling instant, say tk {k >n), we can condense the 
preceding expressions into two sets of equations involving the functions fi,g'^'' 
and hrq as follows: 

n 

J2 fMtk-i - tj) + gl'Lj = {j = k,...,k-n + l) (21) 

i=0 

n 

J2 fihrg{tk-i - tj) = {j = k-n,...,0) (22) 

i=0 

(r = 1, . . . ,p), {q = l,...,m). 

Note that, at time tk, the functions fi, g^'^ will depend on the sampling 
interval lengths {Tk-n+i, ■ ■ ■ , Tk) and so on. 

Now, we multiply (21) and (22) by (wp {q = 1, . . . ,m),{j = k,k - 1, . . . ,0), 
respectively, (uj being the qth impulse input of the system at the sampling 
instant tj) and summing all these expressions, we get 
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m k — 1 

k—n n—1 

+ /n ( E " ^')"?) + E (23) 

1^0 3=0 

m k 

= ^^Kq(tk - U)u1. 
q=l 1=0 

Then, making use of the convohition expression 

m k 

yl = Y.ll^r,{tk-ti)u'> (24) 

q=l 1=0 

{yl being the rth output of the system at time tk) the above expression 
becomes 

n m n — 1 

yi-Y.f^yi-^ + llll97<-, (25) 

1=1 q=l j=0 

which is the input output description for hnear time-invariant MIMO sys- 
tems sampled in an aperiodic way. Each system output at time can be 
written as a hnear combination of the same output and of the different inputs 
at previous instants. The expression (25) generahzes to the aperiodic case the 
well-known input/output representation for linear systems sampled periodically. 
The sequence of sampling intervals is implicit in the arguments of the functions 
/i, (/J*. Consequently, this freedom in the choice of the sampling instants can 
be used in the solution of control problems, propagation of measuring errors, 
parameter estimation, and related topics |11|-|14|. 

It is convenient to note that the functions fi are the same for every system 
output, while the functions gj'' depend on the corresponding impulse response 

3.2 Simplified Computation of the Functions fi, 

Companion matrices are an important example of nonderogatory matrices, which 
have only one (normalized) eigenvector associated with each distinct eigenvalue. 
This means that 

(i) the Jordan canonical form is clearly simplified (there is only one Jordan 

block for each distinct eigenvalue); 

(ii) the similarity transformation of the given matrix to the Jordan form can 
be obtained in a standard way. 



6 



Thus, the computation of the Jordan canonical form for this kind of matrix 
is quite easy. Indeed, 

A ^ TJT^i (26) 

where J is the Jordan canonical form of the matrix A and T is an invertible 
matrix of a well-known general form |8.^. In this way, (14) becomes 



1=0 

with 



'^fn-i{z)xQexp{Jai) (27) 

1=0 

ai = zi + ... + zi (Z = l,...,n),(ao = 0) (28) 



xo = CT. (29) 

Let (ipi) {I ~ l,...,n) be the fundamental system of solutions of a nth 
order homogeneous linear differential equation whose characteristic polynomial 
is d{s). 

In this case, factorizing Det{xo exp{Jai)) and cancelling common factors in 
(15), the general form of the functions fi can be simplified to 

f,^^ {t = l,...,n) (30) 

with 



A = Det{ipi{a,)) (.? = 0, . . . , n - 1), (? - Z, . . . , n) (31) 

and Ai is analogous but replacing the argument of the ith column by a„. 
The functions g^'' can be easily computed from (21) for the new simplified form 
of the functions fi . 



3.3 Choice of the SampHng Sequence 

The procedure developed imposes nonrestrictive conditions on the sampling 
sequence in order to guarantee the linear independence of the vectors {Go{z), . . . , 
G„_i(^)). 

Strategies to determine the set /" C whose elements z verify the above 
condition can be found in 113J (in an analytic way) and in |14| (in a geometric 
way). 

4 An Illustrative Example 

Let H(s) be a two-input, two-output transfer function matrix 

(32) 

(s+l)(s+2) 7+2 I 
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d{s) = s^+3s + 2 = {s + l){s + 2) in = 2) 



(33) 



(pi{t) = exp{-t) 



(34) 



ip2{t) = exp[-2t) 
The impulse response matrix will be 



m = 



(35) 
(36) 



exp{—t) 2exp{—t) 
'exp{—t) + exp{—2t) exp{—2t) 

In this case, for an arbitrary instant tk we know that [13) tk G (^fe-i, oo) and 
so on. 

Thus, choosing Tk-i — 0.8, Tk — 1.1, the expression (25) can be computed 
as follows. 

According to (31) A = 0.24743, Ai = 0.12719, As = -0.02014. 
Consequently, /i = 0.51407 and /2 = -0.08142. According to (10), the 
functions g^"^ will be 



g'o' g'o' g? 

?1 ?1 

5o 5i 5o 5i 



1. -0.18121 2. -0.36241 
0. -0.22207 1. -0.40327 



(37) 



Substituting the previous expressions, we write the input/output relations. 
,1 



y\ ] = 0.51407 ( 

Vk J \ Vk-l 



0.08142 ( ^^-2 I + 

yk-2 



-0.18121 2. 
-0.22207 1. 



-0.36241 
-0.40327 



^,2 



V 



(38) 



For each new sampling instant, the functions /i, g^'^ must be computed 
again. The lengths of the sampling intervals can be chosen in order to optimize 
a particular performance criterion. Further, difficulties may arise in the practi- 
cal implementation of equidistant sampling as, e.g., the equidistance might be 
disturbed. The formulation developed above can be used to pursue the propa- 
gation and consequences of this inexactitude. 

Finally, it should be mentioned that the use of well-known numerical meth- 
ods for the problem of the optimization of aperiodic sampling instants leads to 
good results in concrete cases as it may be seen in [11]- [13]. The formulation 
proposed allows us to use these methods for an I/O modeling technique. 
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5 Conclusions 



The external description developed provides a system model well adapted to 
the real experimentation conditions although presents the limitations inherent 
to the use of the transfer function. 

The formulation considered emphasizes the importance of the sampling se- 
quence against other system parameters. 

The particularization to the periodic case is immediate and represents an 
alternative to the classic discretization methods without using the Z-transform. 

The proccckue is simple and no use of polynomial matrix theory is required. 

This I/O modeling technique allows us to choose the sampling instants in 
order to improve the numerical aspects in problems such as identification, con- 
trol, propagation of measuring errors, . . ., etc. Its use is merely a question of 
an appropriate choice of the performance criterion. 
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